Drag reduction by polymer additives from turbulent spectra 
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We extend the analysis of the friction factor for turbulent pipe flow reported by G. Gioia and P. 
O ' Chakraborty (G. Gioia and P. Chakraborty, Phys. Rev. Lett. 96, 044502 (2006)) to the case where 

vN ^ drag is reduced by polymer additives. 
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INTRODUCTION 



r^ ' When a fluid flows through a pipe of circular section and radius R it experiences a pressure drop per unit 

^~^, length of pipe dp/dx — — 2ro/i?, where tq is the stress at the wall, tq has units of energy density and is 

'"O ' commonly parameterized as in the Darcy-Wcisbach formula 

i 



ro = ^py (1) 



where p is the density of the fluid and V the mean velocity. The coefficient / in eq. ([T]) is the so-called 
friction factor [lH5|- 
r^ ' For a given pipe, the friction factor is a function of Reynolds number Re 

&.: 

en: R«^^ (2) 

> , 

^D , where i^ is the kinematic viscosity of the fluid (as distinct from the dynamic viscosity /i = piy) . It presents 

^^ ' three power law-like regimes separated by transition regions. For laminar flows (Re < 10"^), / = 64/Re; for 

— developed turbulent flows (10'^ < Re < 10^) it obeys the Blasius Law / = 0.3164/Re ' , and for larger values 

of Reynolds number it converges to an asymptotic value determined by the pipe roughness. 
|0 I G. Gioia and P. Chakraborty [6[ have presented a theoretical model whereby these features of the friction 

arc easily derived from the Kolmogorov spectrum of homogeneous, isotropic turbulence. See also |7H14|. We 
shall refer to this as the momentum transfer model (MTM) 

An immediate prediction of the MTM is that, in situations where the turbulent spectrum deviates from 
the Kolmogorov form, the friction factor should also deviate from the Blasius Law. This prediction has been 
confirmed in the analysis of two dimensional flows |ll| . In this note, we wish to perform a similar analysis 
}J] ' for three dimensional pipe flow in the presence of drag reducing polymer additives. 

5^ I It is well known that adding a few parts per million of certain polymer additives to a fluid causes a drastic 

reduction in the friction factor |15l-l26| . Use of this effect in improving the efficiency of oil and natural gas 
pipelines is widespread. There are indications that the phenomenon is not confined to pipe flow, but that 
the presence of the additives affect turbulence even in the homogeneous and isotropic limit [27[ . 

In this paper we will show that, given a turbulent spectrum for the pure solvent consistent with both the 
Blasius Law / = .31/Re^' and the Poiseuille Law / = 64/Re, then there is a deformation of this spectrum 
that reproduces the phenomenology of drag reduction, both in the asymptotic universal limits and with 
respect to the concentration dependence. 

To incorporate the polymer to our model we shall adopt the theoretical framework provided by the 
so-called finitely extensive nonlinear elastic model supplemented by the Peterlin approximation (FENE-P 
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model) [28|,|29[. To obtain the turbulent spectrum under homogeneous isotropic conditions, we shall map the 
nonlinear equations of the FENE-P model into an equivalent stochastic linear system, constructed to produce 
the right turbulent spectrum in the zero polymer concentration limit. By solving this linear problem, we shall 
find the spectrum as modified by finite polymer concentration, and we shall show that indeed it becomes 
universal in the large concentration, large Reynolds number limit. Moreover, by adopting the MTM, we will 



derive a power law dependence for the friction factor / sa Re^^", close to the experimental result 15] 

This paper is organized as follows: in next Section we review the MTM of the friction factor in the 
absence of the polymer. The original presentation of the MTM made contact with the Blasius and Strickler 
(for rough pipes) asyniptotic regimes, but did not discuss in any detail the transition from the Blasius to 
the Poiseuille regimes [g, [ij] . To incorporate the Poiseuille regime within the MTM framework, we analyze 
the flow into a central region, where velocity fluctuations play an important role, and an outer region where 
fluctuations are negligible. We apply the MTM prescription to find the Reynolds stress at the boundary 
between these two regions, and then relate it to the stress at the wall by solving the Navier-Stokes equation 
in the outer region. We show that the resulting model gives the Poiseuille Law at low Reynolds numbers. 
At high Reynolds numbers the flow in the central region may be described as a Kolmogorov cascade, and in 
this case the model yields the Blasius Law. 

In the following Section, we introduce the polymer. We show in the Appendix that the effect of the polymer 
in the outer region is negligible. To find the flow in the central region, we approximate the Navier-Stokes 
and FENE-P equations by a linear system driven by a stochastic force. The stochastic equation for the fluid 
is chosen by requiring that, in the absence of the polymer, it reproduces the spectrum of fluctuations as 
described in refs. 0, [iJl- We then add the coupling to the polymer stress tensor as dictated by the FENE-P 
model. The equation for the polymer deformation tensor is a Hartrec approximation to the original FENE-P 
equation. 

The model leaves several parameters indetermined, the most important being the relaxation time of the 
polymer. We determine this parameter by requiring that at high Reynolds number the relaxation time for 
the polymer is proportional to the revolving time for the eddies that dominate momentum transfer in the 
MTM. These are the eddies whose size matches the width of the outer region. At low Reynolds numbers, 
the relaxation time regresses to its equilibrium value. Since the model is not sensitive to the details of the 
relaxation time dependence with respect to Reynolds number, we assume a simple interpolation formula that 
yields the proper asymptotic values. 

The result of this analysis is a friction factor - Reynolds number dependence containing five dimensionless 
parameters. We determine these parameters matching to experimental results, namely the Poiseuille and 
Blasius laws without the polymer, the Virk asymptote [l5| at large polymer concentrations, and finally the 
detailed data presented in [l5| for finite concentration. Having obtained a suitable set of parameters, we 
display the results in Section IV. 

We conclude with a few brief remarks. In the appendix we discuss the FENE-P model in the outer region. 

2. THE MTM APPLIED TO THE PURE SOLVENT FRICTION FACTOR 

In the absence of the polymer, the dynamics of the solvent is described by the Navicr-Stokcs equation 

d 



s-'"^' 



UP + Vq [upu«] + ypp = (3) 



We arc interested in a stationary flow within a straight pipe of circular section and radius R. Let z be 
the coordinate along the pipe. The flow may be decomposed into mean flow and fluctuations as U^ = 
U (r) zP + uP, where z^ is the unit vector in the z direction. We call V the average value of U across the 
section of the pipe. We define a dimensionless radial coordinate ^ = r/R. 

We are interested in a high Reynolds number regime where the mean flow is very flat in the central region 
of the pipe, from ^ = to ^ = ^* = 1 — 6*, say, and there is an outer layer from there up to ^ = 1. The 
precise form of the velocity proflle in the central region is not a critical concern and we shall take it as simply 
flat, with amplitude U* . We adopt the convention of computing the Reynolds number as if the central flow 



filled the whole pipe; this is only a matter of convenience. The central flow is characterized by a Reynolds 
number 



Re* 



2U*R 



In the outer layer we neglect the fluctuating velocity u^. Then the Navier-Stokes equations reduce to 



d 1 d dU _ 



The solution that vanishes at ^ = 1 reads 



(4) 



(5) 



U[i]^U*[a{l-e)+bln[ 
Asking the mean velocity profile to be continuous we get 



(6) 



1 ^ a {1 - C^) + blii[C 
We also have the shear stress at ^*, namely 



(7) 



where 



Tturb = r"'^ 



1 dU ,^^, „ ^, b 

r = 2ar — 

U* dC C 



(8) 



(9) 



We can write the constants a and b in terms of x* a-s 



A[e 



b = 



r 



[i+rinrixi 
[2r-(i-r')x*] 



A[e 
A[r] = 1 - r^ + 2r^ In [r] 

Our interest is to find the average velocity, which enters in the Reynolds number eq. ([2]) 



V = C^U* + f dxU [yM = U*h [6*] 



(10) 



(11) 



h [S*] = 



2A[e] 1 



+ x £. 



i + (i+r 



ln[e 



(l_^*2 



and the stress at the wall 



(12) 



pvU* 



(13) 



where 



1 dU , , 
X--j^^[l]=2a-b (14) 

If we parameterize tq as in eq. ([T]), then the friction factor 

8to _ 64 xU* 
■' " pV^ ~ Re W ^^^' 

If ^* — > then a— >1,6— !>0,V^— >[/*/2 and x — > 2, so we recover the Poiseuille Law. 

In the general case, we need to relate \J* and x* to ^* to obtain the friction factor - Reynolds number 
dependence in parametric form. These relations are provided by the momentum transfer model (MTM). In 
the central region, each scale 5 is associated to a velocity u [(5, [/*] . The MTM claims that 

a) the width 8* of the outer layer is also the Kolmogorov scale of the central region, namely the scale at 
which the Reynolds number is 1 

^^l^iM = 1 (16) 

Note that in the original presentation of the MTM only a proportionality between 8* and the Kolmogorov 
scale is required p . We have adopted the more restrictive criterion (fT6)) to simplify the discussion below. 
Let us write 

u\6\U*\^U*v\8*\ (17) 

then cq. p6p may be rewritten as 

Re* = j:^ (18) 

Q*v \a*\ 

b) the shear stress at the boundary of the central region is 

Tturh^pU*uW,V*\ (19) 

As in case (a), we have opted for postulating an equality where the original MTM only asks for proportionality 
S. Therefore 



X* = ^ (20) 

and the constants read 

1 r f* 

^-^[i-r] (21) 



leading to 



We may now write the friction factor - Reynolds number dependence in parametric form 

Re = C[S*] 
f = V[S*] (23) 

where 



5*v [(5*1 



V>m = ^^W<51 (24) 

h [o*\ 

In this formulae we already have explicit expressions for h [S*] (cfr. eq. ([T2|) ') and xo [S*] (cfr. eq. p2|)'). but 
we need a detailed model of the velocity fluctuations in the central region to derive v [S*]. This shall be our 
concern in the rest of the paper. 

Wc have already remarked, however, that if 5* — > 1 then the parametric equations (j23p reproduce the 
Poiscuille law, and it can be seen by inspection that a Kolmogorov scaling v [S*] oc S*^^^ when i5* — >■ 
will produce the Blasius law, if appropriate values for the several constants in the theory may be found. 
Therefore we may be confident that our model successfully reproduces the limiting behaviors. 

3. SOLVENT - POLYMER INTERACTION 

3.1. The FENE-P model 

In this section we consider the modifications of the above picture due to the addition of the polymer. We 
shall adopt the so-called FENE-P model. The fluid velocity obeys the incompressibility condition VpU^ ~ 
and a modified Navicr-Stokes equation 



9 ^2 



\JP + Vq[\5PW] + VPP^-Vq{ppTP'i) (25) 



where P is the pressure, p the fluid density, pp the polymer density and T^'^ a polymer stress. T^' is modeled 
in terms of the polymer deformation tensor A^' as 

TP'' = u:),,,h{K)K^'^ (26) 

where w/ree is the frequency of free oscillations of the molecule, A = (l/3)trAP'' and h is some function that 
is close to one under equilibrium conditions A ~ K^q and diverges as A approaches maximum elongation 

The evolution of the deformation tensor is determined by the drag from the fluid and the polymer elasticity. 
Neglecting the inertia of the molecule, we get 



l + U^^'- 



AP9 = (V^U^) A''« + A^'' (V^U«) - ts {Ti"i - TPJI) (27) 



ts is the time scale in which a freely moving bead from the polymer would come to rest with respect to the 
fluid. 



3.2. Flow in the central region 

We shall make the approximation that the flow in the outer region is not affected by the polymer. This 
issue is further discussed in the Appendix. Under this approximation, the analysis in the previous section 
remains valid, and the only effect of the polymer is changing the functional form of u [6, U*] in eqs. ()16|) and 
([TO)) . To find this, we only need to consider the fluctuating part of the velocity. We shall consider only the 
homogeneous, isotropic case, since the behavior of the fluid in this case determines the friction factor in the 
MTM. To take advantage of the symmetries of the problem, wc shall decompose the deformation tensor into 
its scalar and traceless parts 

AP9 ^ A(5P« + XP" (28) 

A^ = 0. Moreover, we shall assume that A is both space and time independent. The stress tensor is 
decomposed into a similar way 

TP1 ^ TS^" + uj}^^^h (A) AP« (29) 

where 

T = u%,,h[K)k (30) 



Taking the trace of eq. (p7|) we obtain 

3ts [T - Teq) = 2AP'VpU, (31) 

Subtracting the trace from eq. (|27)) we get 



dt 
where 



O 1 

AP« = A (VPu« + V«uP) - -AP9 + SP'i (32) 



r 



1 .2 .>.^,__Tts 



T 



w^_/i(A)i5 = -p (33) 



A 



and 



SP" = (V,.uP) A""' + A'^'' (V^u«) - ^A"^" (V^u,) (5^9 - u'-V^A^' (34) 

whose ensemble average must be zero from the symmetries of the problem. Wc shall neglect 5''"' in what 
follows. 



3.3. Equivalent linear stochastic model 

It is clear that the full Navier-Stokes is too complex for analysis, unless numerically [20|. To make progress, 
wc shall substitute the Navier-Stokes equation by a linear stochastic one, devised to give the right spectrum 
in the absence of the polymer. 

Let us begin by Fourier decomposing the fluid velocity 



We postulate for the Fourier components a dynamic equation 



"''^•"^/^(^'""""""'^i^> ''^' 



[-IUO + <7fc] < = F^ M (36) 



where F^ is a Gaussian random source with self correlation 



{Fl M F^, [u']) = (2^)*^ 5 (k + k') 5 {uj + uj') Al'N, (37) 



Ar = '5^^-^ (38) 

A representation like this may be derived from the functional approach to turbulence, where the left hand 
side of eq. (|36|) is identified as the inverse retarded propagator, and the self-correlation eq. ((37|) is given 
by a self-energy [2l[ [sol- We shall be content to propose simple expressions for o-fe and Nj, to reproduce the 
known turbulent spectrum. 

In the inertial range, we expect cjk and N^^ to depend on the only dimensionful parameter e, which is the 
energy flux feeding the Richardson cascade. On dimensional grounds |3l| 

ak^iyo{k\y^' (39) 

where vq is a dimensionless constant to be determined presently. The turbulent spectrum Eq [k] (where 
the subscript denotes that this is the spectrum in the absence of the polymer) is defined from the mode 
decomposition of the turbulent energy 



/•OO 

(u^ [x,t]) =2 / dkEo[k] 
Jo 



(40) 



Explicitly 



Eo{k]=e^ (41) 

so we recover the Kolmogorov spectrum Eq [k] = C/^e^/'^fc"^/'^, where Ck ~ 1-5 is the so-called Kolmogorov 
constant [21|, provided 

Nk = ''oCk^ (42) 

The ansatz eq. ([M]) for ak is equivalent to a scale dependent viscosity ly^ = k^^ak = i^oe^^^k^^^'^. Under 
Kolmogorov scaling the velocity associated to a scale S = k^^ is Uk = \/3Ck {^k^^^ . We find the Reynolds 
number of the effective linear theory as RCe// = ui^/kiy^. = ^JZCkIvq. Identifying this with the physical 
Reynolds number of the central region we get 

-0 = ^^ (43) 



This simple picture must be modified to account for the dissipative range. In the dissipative range fluctuations 
are strongly suppressed 



N,.^ 



V^: 



K 



5'Rk 



Re* fc3 



(44) 



with /? Si 1/2 a dimensionless number |14l |. The strong suppression of fluctuations dispenses with further 
discussion of a^ in this range; on general grounds we expect it will approach its bare value t^fc^, but for 
simplicity we shall use the incrtial form eq. ()39|) in the calculations below. 

At very long wavelengths the spectrum must turn over and approach the von Karman spectrum E [k] ex fc^ 
[32l |. To obtain this we add one further factor, turning the noise correlation into 



Nk^ 



v/3^ 



^^17/3^8/3 



^-135'Rk 



Rc* 



-f+iRky 



17/6 



(45) 



This yields the same spectrum as assumed in ref. 



3.4. Solving the effective linear model 

Adopting the linear effective model as a suitable description of the eddy dynamics, and for a constant 
polymer density pp ~ cp, we get, instead of eqs. ((25|) and ([27| . the linear system 



--+-fc]<-^Mr-^n-] 



(46) 



T 



Al^''=*A(kX + kH) 



(47) 



Eliminating A^'^ we obtain a second order equation for u^ (we also use the incompressibility constraint 

KK = 0) 



P[io]ul = 



—iuj + 



KM 



(48) 



r ii 


1 


CTfe + - 





T 


r 



P [lS\ = lo" -\- iuJ 

= (uj — S_|_) (w — s_) 
s± are the free frequencies of the system and are given by 



C7k + — Afc' 



(49) 



1 



-i ( o-fc + - I ± 



/4cAfc2 

Tts 



CTfc — 
r 



(50) 



We see that there are two different flow regimes. When the discriminant in cq. (|50)) is negative, both 
eigenfrequencies are pure imaginary. The imaginary part of both is always negative, so the flow is always 



stable. We shall call this the overdamped regime. This regime prevails in the energy range (where fc — >■ 0, 
and therefore also Cfc) and in the dissipative range, where a^ is very large. 

On the other hand, precisely because Cfc goes from at fc = to a very large value in the dissipative 
range, there must be some interval where a^ ~ 1/t and the discriminant is positive. In this regime the 
free frequencies have nonzero real parts, although they still describe damped oscillations. Wc shall call this 
the undcrdamped regime. As the concentration c grows, the undcrdamped range expands and essentially 
becomes identical with the inertial range. 

From the solution to eq. pS)) and the noise self-correlation eq. ([57)) we identify the spectrum in the 
presence of the polymer as 



E[k]=k'-Nk[^J[k] + J2[k] 



(51) 



where 



and 



Actually these integrals are related 



J[k] 



J2[k] 



dw 1 



TT P[UJ]P[-UJ] 



du) uj 

~ P [uj] P [-uj\ 



(52) 



(53) 



J2 \k\ = 



o-fc 



cA/c^ 



ts 



J[k] 



(54) 



Evaluating the integral 



J[k] 



s+s- {s+ + s_) 



(7k 



[1 + CTkT] 



ts^k 



(55) 



and so the spectrum is 



E[k] 



1 + CTkT - 



"""^ [1+^feT] 



crAfc^ 



ts 



1 + 



tsck 



(56) 



3.5. Identifying the free parameters 

To give meaning to eq. (|56p wc must know the way parameters such as e, r and A depend on S*. The 
determination of these parameters is the subject of this section. 

Let us begin with the expression of u [S*, U*] in terms of the spectrum (cfr. eq. (j40|) ) 



,[S*,U*]^ =2 / dkE[k 



(57) 
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writing x = 5* Rk this becomes 



[5*^ U*f = 2CkS*'^^^ {Ref^^ / dx 



„4„-x/2 



1 + a^T 



rAfc^ 



ts 



[7<5*2+:e2]1V6 j^ 



+ CTfeT 






(58) 



To obtain e, we observe that energy is fed into the Richardson cascade at a scale i?^*. Therefore we expect 



/ \ -i TT* 



with K a dimensionless parameter to be determined. This leads to 



(59) 



K V X 



2/3 



CTfc 



2 i?2 J*2/3^*l/3 



(60) 



To find T, we expect that at large Reynolds numbers r will be proportional to the revolving time for eddies 
of size 5*, namely 



(61) 



u[5*,U* 

where a' is a new free parameter. For lower Reynolds numbers, we expect r will regress to its equilibrium 

A t 

free 



value l/w^^ggtg. To interpolate between these regimes, we assume 



1 



W/ree^S ( ^ + j:^ 



where 



^' a'R^u%^^ts 



A is related to r through eq. ([55]) . If we neglect the ratio t^eql l^max and parameterize 



/i(A) 



then 



A = 



A. 



"5^ 



This leads to 



(62) 



(63) 



(64) 



(65) 



" "O 2/3 

CTkT = a— TTTTo :r-;X ' 



ei/3 51 + (5* 



(66) 



11 



crAfc^ 
ts 



"o 



CO [(52 +5*2]^ 



(67) 



1 ce 



■1/3 



S'o 



cAk^ 



ts(Jk a CO <5*4/3 [^2 ^ g* 



.4/3 



where 



(68) 



(69) 



Co 



From eqs. 



and ([T7)) we get 






(70) 



v[5*] 



J*V3 




dx 



4 ■ 

a; e 



-x/2 






„2/3 



+ 



CO [52+5-2 



|T2; 



1/2 



[7(5*2 



.2117/6 r 



1 



C-i/3 52+5*2 



,2/3 



1 C 4*1/3 52 

a CO 5W3 ^52 .j.^. 



,4/3 



(71) 



We see that the solution depends on the parameters n, 7, a, cq and 6q. In principle, each of these could be a 
function of Reynolds number or other dimensionless combinations, This would turn the parametric relations 
eq. (j23p into implicit equations. For simplicity, we shall model them as constants. 

Let us check that our model shields the appropriate limiting behavior. To begin with, to obtain Re ^ 
when (5* — ?> 1 we need that v [S*] should diverge in this limit, which indeed it is a result of eq. (|7ip . This does 
not contradict the fact that U* remains finite because in this limit the scale R6* is larger than the scale R^* 
at which energy is injected into the fluctuations. In this limit, of course, the fluctuations themselves cannot 
be regarded as turbulent in any conventional sense, and our model should be regarded as an extrapolation 
which "saves the appearances" . 

If c = and (5* -^ then eq. dTT]) predicts v [S*] ex 6*^^^ and leads to the Blasius Law. We see that 7 7^ 
is necessary to obtain the "hump" feature in the friction factor plot 0, [ij] . 

If (5* — > but c is large, then eq. (|7ip predicts v [S*] ex S* and then the parametric relations become 
/ oc Re~^' ^, in reasonable agreement with Virk's asymptote p^ l33l - [36| . We see however that even for large 
concentrations drag reduction will be very small as long as S* 3> ^0 



RESULTS 



In this section, we shall use the previous analysis to obtain concrete estimates of the friction factor. We 
adopt the following values for the free parameters: k = 0.02, 7 = 40, a = 0.05, 1/Sq = 5000 {60 = 0.014) 
and Co = 43.6/0.00035 = 124600 wppm 

In fig. ([T]) we show the concentration dependent spectra for Re :=: 10^ and c/cq = 0, 0.01 and 1. 

We now discuss the solution to the parametric equations (P5)) . In fig. ^ we show the friction factor for 
the pure fluid. The transition from the Blasius to the Poiseuille regimes is clearly seen. To obtain a more 
accurate fit to experimental data would require the introduction of a more complex spectrum and is not 
relevant to the discussion of drag reduction. 
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FIG. 1: [Color online] The concentration dependent spectra for Re — W . The full lines, from the top down, 
correspond to c/co = 0, 0.01 and 1; the concentration scale cq = 124600 wppm is defined in the text. The dashed 
line represents the Kolmogorov spectrum. We also show the beginning and the end of the inertial range in the pure 
fluid limit 
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FIG. 2: [Color online] The friction factor in the absence of the polymer (full line). The straight dashed lines represent 
the Blasius (short dashes) and Poiseuille Laws (long dashes) 



In fig. dS]) wc add, to the c/cq = line in fig. ^, the friction factor dependence for c/cq = 0.01 and 1. 
We also add the Virk asymptote for comparison. 

It is convenient to introduce the Prandtl-von Karman variables X = In [Re-y//] ^^d Y = l/\/J- In these 
coordinates the Poiseuille Law becomes Y = 6^^/64. In the turbulent regime, Y is described by the Prandtl 
Law Y — 0.81X — 0.8. While the Prandtl and Blasius Laws, which reads Y — (.31) e^^'^ , have very 
different mathematical expression, they give equivalent results in the range of Reynolds numbers we are 
considering. 

In fig. 131) we reproduce fi g.(|3l) in Prandtl-von Karman coordinates. We have also added the Virk asymptote 
Y = 4.12X - 19.06 [ll,[33,l3l. 

In fig. ^ wc compare the result from our model to the experimental data presented in ref. |17| 

A comment is in order about the value of cq used to draw these plots. It is clear that the value of cq we 
are using corresponds to a very high concentration, probably higher than any used in actual experiments 
[l9|. However, our cq does not represent the onset concentration. As shown in fig. ^ and will be seen again 
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FIG. 3: [Color online] The friction factor for non-zero concentration. The full lines, from the top down, correspond 
to c/co — 0, 0.01 and 1. The dashed lines, from the top down, correspond to the Blasius, Virk and Poiseuille Laws. 




X 



FIG. 4: [Color online] The friction factor from Fig. eq. ((3]) in Prandtl-von Karman coordinates. The dashed lines 
are, from the bottom up, the Prandtl, Virk and Poiseuille Laws 



in the following figures, substantial drag reduction is seen at a concentration of 10^'^ cq, and so a high value 
for Co is to be expected. This said, it is clear that the multiplicity of parameters and the lack of independent 
derivation and/or determination of at least some of them is a weakness of our model and an area for further 
work. 

To analyze the dependence of the friction factor on concentration, it is convenient to introduce the frac- 
tional drag reduction Rp 
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FIG. 5: [Color online] The friction factor for small polymer concentration. The symbols represent the data presented 
in Fig. 2a of ref. [17| . They correspond to measurements of the friction factor in solutions of a single polymer (N750) 
in five different concentrations: pure solvent (full circles), and 43.6 (triangles), 98.6 (squares), 296 (circles) and 939 
(diamonds) w. p. p.m. .For the theoretical lines we have used the parameters k — 0.02, 7 — 40, a = 0.05, 1/5^ = 5000 
{So = 0.014) and cq = 43.6/0.00035 = 124600 wppm. Therefore the lines in the plot correspond to c/cq = 0.00035, 
0.00079, 0.0023 and 0.0075 



Rf = 1 



/o 



(72) 



Rf ^- as c ^ by definition. When c — ?► 00, on the other hand, it reaches a finite asymptotic value 
We also introduce the intrinsic drag reduction Ri 



Ri 



Ri 



(73) 



Rj goes to zero when c — > 00, but when c — > it reaches a finite value Rifi- Following Virk [I3I, we define 
G = Rific/ RF,max and D = Ri/Ri^. Then the following empirical relation holds 



D = 



1 



1 + G 



(74) 



We plot these quantities in figs. (jS]), (O and ([5]). The agreement of this last figure to eq. ([71)) is remarkable. 
This shows that our model not only predicts there will be a maximum drag asymptote, but also reproduces 
the full concentration dependence. 



FINAL REMARKS 



In this paper we have shown that, given a suitable expression for the turbulent spectrum in the absence 
of the polymer, there is a deformation of it that reproduces both the maximum drag reduction asymptote 
and the concentration dependence of the intrinsic drag reduction. 

Our treatment is admittedly not a self-contained derivation of the friction factor; for once, the model 
allows for a large number of parameters which are not determined independently, but simply chosen to fit 
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FIG. 6: [Color online] The fractional drag reduction for Re = 10"" 
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10"" 1 ^0 

FIG. 7: [Color online] The intrinsic drag reduction (multiplied by co) for Re = 10^. 
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10"" 1 10 
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FIG. 8: [Color online] D vs G for Re = 10^. The full curve is derived from our model, the dashed curve is the 
empirical relationship eq. (|74ll . 
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experimental data regarding the friction factor itself. Granted this, we believe each step in our argument 
is well motivated. The discussion of the polymer, for example, is based in what essentially is a Hartree 
approximation to the FENE-P model, and as such stands on a well trodden (theoretical) path. 

One of the most remarkable features of the drag reduction phenomenon is the universality of the maximum 
drag reduction asymptote. In our model, universality obtains from the fact that we assume that, when 
Reynolds numbers are high, the polymer relaxation r is determined by the eddy revolving time alone, 
independently of polymer characteristics (which do play a role at lower Reynolds number). This assumption 
is inspired in Lumley's "time" criterion [16|. However the Lumley criterion refers to the characteristic time 
pv /tq, and also the way the polymer characteristic time is defined is different from ours. 

The quantitative fit to the Virk asymptote depends on the other hand on the assumption that the lifetime 
of a velocity fluctuation grows with Reynolds number as in eq. (PH)) . with the dimensionless factor eq. (|43p . 
This condition follows from the requirement that both the physical and the equivalent (linear) flows share 
the same Reynolds number at large scales. Once it is accepted, it follows that the random driving must also 
weaken with Reynolds number, as could be expected from fluctuation-dissipation considerations [37[. 

Overall, we believe the results of this paper are a success for the MTM, complementing earlier studies of 
the friction factor in two-dimensional turbulence [ll[. We offer them as a simple theoretical template for 
more fundamental approaches. 
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Appendix: The FENE-P model in the outer region 

In this appendix we shall discuss the arguments behind the contention that the polymer does not affect 
the flow in the outer region. We thus neglect the fluctuating velocity and assume U^ = U (r) z^, P = ~Pqz, 
and A^'' = A^'' (r) . The left hand side of eq. (|27)) vanishes and we get six algebraic equations 

ts (T^^ - T^^') = 2C/'A^'- (75) 

where U' = dU/dr 

tsT''' = t/'A'"'' (76) 

tsT""^ = U'A"'' (77) 

Therefore 



writing 



^0r = A"^ = (79) 



TP" = oo%,,K,,9^'^ (80) 
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where g^"^ = diag (l, 1, l/r^) we get 



A"' 


h{A) 


, ee 


Aeg 




r2/i (A) 




U'Aeg 



(81) 
(82) 
(83) 



A- = -11^ + ''^^ , (84) 



J 2U'^A 



Plus the consistency condition 



In the limit when A w A„ja^ ::!> Ae„ eq. ([SS]) reduces to 



/i (A) 1 = ., / 2 (85) 



,3 2U'^Aeg 



h (A)^ = , , "l (86) 



Given this form of the deformation tensor, the only nontrivial Navier-Stokes equations is the z equation, 
which yields 

dS,S.dCd£, 
where 

^ = ^ + 2^T- (88) 

dt, dt, V 

We sec that this is the same equation as without the polymer, only the fluid velocity gradient is "corrected" 
by a term 

^-^^/.eefe(A) ^^,,^c A,, dU ^ C Ae, dU 

V V tgh (A) di aco A„,axh (A) d£, 

Even if 1/a = 20, in the relevant regime all c/cq, A^q/ A^ax and 1/h are very small. So the correction to 
dU/d^ may be disregarded for all practical purposes. 
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